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GAUSSIAN LIMITS FOR RANDOM MEASURES IN GEOMETRIC 

PROBABILITY 

By Yu. Baryshnikov and J. E. Yukich 1 
Bell Laboratories and Lehigh University 

We establish Gaussian limits for general measures induced by 
binomial and Poisson point processes in d-dimensional space. The 
limiting Gaussian field has a covariance functional which depends 
on the density of the point process. The general results are used 
to deduce central limit theorems for measures induced by random 
graphs (nearest neighbor, Voronoi and sphere of influence graph), 
random sequential packing models (ballistic deposition and spatial 
birth-growth models) and statistics of germ-grain models. 


1. Introduction. The purpose of this paper is to provide a methodology 
for showing that renormalized random point measures in geometric proba¬ 
bility converge weakly to a generalized Gaussian field. We focus on random 
point measures, defined on the Borel subsets of R rf , of the following types: 

(i) point measures associated with random graphs in computational ge¬ 
ometry, including nearest neighbor graphs, Voronoi graphs and sphere of 
influence graphs, 

(ii) point measures arising in random sequential packing models, includ¬ 
ing random sequential adsorption (RSA) and spatial birth-growth models, 
and 

(iii) point measures associated with germ-grain models. 

The total mass of random point measures yields random functionals, 
which in the context of the measures (i)-(iii), have been extensively studied; 
see [2, 21, 23, 24, 26, 32], [5, 7, 8, 9, 11, 15, 25, 27] and [12, 13, 22, 26], re¬ 
spectively, as well as the references therein. With the exception of [13], the 
study of the random measures (i)-(iii) has received considerably less atten¬ 
tion. We show here after renormalization that measures of the type (i)-(iii) 
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converge to a generalized Gaussian field; that is, their finite-dimensional 
distributions, as described by the action of the measure on continuous test 
functions, converge to those of a generalized finitely additive Gaussian field. 
The results relate the large-scale Gaussian limit properties of renormalized 
random point measures to the small-scale properties of the underlying bino¬ 
mial or Poisson point process. 

The general approach taken here, which employs stabilization of func¬ 
tionals and coupling arguments, has the particular benefit of describing the 
limiting variance over large sample sizes as a function of the underlying 
density of points. A similar approach is used in [26], which treats the easier 
problem of finding limiting means. 

Random measures considered here assume the form £(®; X)5 x , where 

A is a random point set in R rf , S x is the Dirac point measure at x and £(x; X) 
is a weight representing the interaction of x with respect to X and is usu¬ 
ally defined in terms of the underlying geometry. For all constants A > 0 and 
probability densities k, let V\ K be a Poisson point process with intensity mea¬ 
sure Xk : M. d —> M + . Define the “binomial” point process X n : = {Xi,..., X n }, 
where Xj, z > 1, are i.i.d. with density k. All of our results follow from general 
central limit theorems (Theorems 2.1, 2.2, 2.5) which show that renormalized 
measures of the type 

(1.1) A" 1 / 2 Y £(A 1/d x;A 1/d V XK )8 x , A > 1, 

as well as their respective renormalized binomial counterparts, 

(1.2) n" 1/2 Y £,(n 1/d Xi-,n 1/d X n )5 Xi , n> 1, 

converge weakly as A —* oo (resp. as n —> oo) to a Gaussian field with a 
covariance functional described in terms of the weight £ and the underlying 
density k of points. 

The general central limit theorem (CLT) for the measures (1.2) implies a 
CLT for the “total mass” functional J2xi£X n £(n 1//rf Xj; v}^ d X n ). k need not 
be uniform and £ need not be translation invariant, showing that even in the 
functional setting, we extend and generalize previous results [2, 3, 5, 24, 25]. 

The proofs are based on the method of cumulants, which requires showing 
that the cumulants of the integrals of the rescaled measures (1.1) against a 
large class of test functions converge to the cumulants of a normal random 
variable. An important tool is “stabilization” of functionals, used heavily in 
[5, 24, 25, 26]. Stabilization guarantees that the pair correlation function for 
the weights £(x,V\), x € M , decays fast enough to prove convergence of the 
cumulant measures associated with (1.1). To show convergence of the first- 
and second-order cumulant measures against test functions, we rely upon the 
“objective method,” which exploits the fact that if £ is locally determined in 
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a sense to be made precise, then the large A behavior of £(A 1//d :r, X^ d V\ K ), x 
fixed, is approximated by the behavior of £ on homogeneous Poisson point 
processes. This idea was developed in [26], a law of large numbers (LLN) 
precursor to the present paper. To show convergence of the higher-order 
cumulant measures, we employ cumulant expansion techniques [20]. 

2. Main results. 

2.1. Terminology. Before stating our main results we introduce some 
terminology similar to that developed in [5, 24, 25, 26]. Let X C be finite 
and y + X := {y + x : x E X} for all y E R d . Given a > 0, let aX := {ax : x E 
X}. For x E I d , |x| denotes the Euclidean norm and B r (x) denotes the 
Euclidean ball centered at x of radius r. u>d denotes the volume of the unit 
ball in M d and 0 denotes the origin of W l . 

Throughout, A denotes the family of compact, convex subsets A C M f; 
with nonempty interior. Let A! denote A together with the space M. d . For 
Ae A f , C{A ) denotes the continuous functions /: A —> R. For / E C(A) and 
fj, a Borel measure on B(A) we let (/,/t) := f A fdu. Given f :M. d —> R, let 
Supp / be the closure of {x E R d : f(x) A 0}. 

Let £(x; X) be a measurable R-valued function defined for all pairs (x, X), 
where x is an element of X. For the moment, we assume that £ is translation 
invariant, that is, £(x; X) = £(x — y; X — y) for all y G W l . When x ^ X, we 
abbreviate notation and write £(x;X) instead of £(x;TUx). 

Any finite X induces the point measure Z]xga £( x 'i X)5 x - For all A > 0, 
let £a(x;T) := £(A 1 / d x; \ l / d X). A density k with support on A G A and a 
weight £ generate scaled random point measures 

lAn : = 

The centered version of A\ K is A\ K := A\ K ~~ ^A\ k i where for all Borel 
sets B C A, E[h\ k (B)] = A j B E[£_>, {x\V\ K )]n(x ) dx. This paper develops a 
methodology for establishing convergence of the finite-dimensional distribu¬ 
tions of the renormalized random point measures A -1 / 2 /if , A > 1. Previous 
work [25, 26] developed laws of large numbers for the total mass functional 

(2-1) /4(A) := E 6, (x-V XK ) 

as well as CLTs [24] for translation-invariant functionals on uniform point 
sets which are “locally determined.” The following concept of stabilization 
makes precise the idea of “locally determined.” For all 0 < a < b < oo, let 
T{a,b) consist of all f :M. d —> M + having support in A! and such that the 
range of / is in [a, b] U {0}. The common probability space (Q,iF,P) for all 
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Vf, f E V(a, b ), can be chosen as the probability space of the Poisson point 
process V* having intensity 1 on R d x R + such that Vf = 7r E d {V* n {(x, h ) E 
U. d :h < f(x)}), where 7r K d denotes projection from R d x R + onto R rf . For 
all r > 0, let V T denote a homogeneous Poisson point process on R d with 
intensity r. 

Definition 2.1. The functional £ is stabilizing if for all A E A', 0 < a < 
b < oo, A > 0, and x E A A, there exists an a.s. finite random variable R(x) := 
R(x,X,a,b,A ) (a radius of stabilization for £ at x) defined on (Q,P,P) such 
that for all / G P(a,b), with Supp/ = XA, and all finite X C XA\B R (x ) we 
have 


£0g (Vf n B r (x)) U X) = £(x; Vf n B R (x)) 

and moreover sup xeR d P[R(x, A, a, b, A) > t] —>• 0 as t —> oo. When £ stabi¬ 
lizes, then for all iGR'* and all r > 0 we define 

f(x\V T ) := lim £(x; V T n B t (x)). 

I —>oo 

Thus R := R(x, A, a, 6, A) is a radius of stabilization if the value of £(x; Vf), 
f G P(a,b), is unaffected by changes outside B R (x ). One might expect that 
exponential decay of the tails of R implies exponential decay of the corre¬ 
lations of £ and thus convergence of A~ 1//2 iz| K ,A > 1, to a Gaussian field. 
This loosely formulated idea figures prominently in interacting particle sys¬ 
tems on the lattice, and also in cluster expansions and the moment method 
in statistical physics [20]. Assuming neither translation invariance of £ nor 
spatial homogeneity of points, we will show that this idea also works well 
in the continuum, where it yields convergence of (/, f €= C(A), 

to a Gaussian field whose covariance depends on the density of points. This 
motivates defining uniform tail probabilities for the radii R(x, X,a,b, A): 

r(t) := r(t,a,b, A) := sup P[R(x, X,a,b, A) >t\. 

x£R d ,\>0 

r(t ) quantifies the region of influence of points in the Poisson point sets Vf, 
whenever / G V(a, b ) and Supp / is a scalar multiple of A. £ is exponentially 
stabilizing if r(t) decays exponentially in t for all a, b and any A G A!. £ is 
polynomially stabilizing if for all a, b and A G A 1 we have (x(t)) 1 ^ 2 t d ~ l dt < oo, 
which readily implies the rough estimate r(t) = o(t~ 2 ). 

The next condition is used frequently in the scaling limit analysis of ran¬ 
dom fields on lattices (e.g., page 193 in [20]) and it is only natural to use 
it in the continuum setting as well. Here and henceforth k is a probability 
density which is continuous on its support and SuppK G A. Let C denote the 
collection of finite point sets in R rf . 
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Definition 2.2. £ has a moment of order p > 0 with respect to n if 

sup K[\^x(x-,V\ k U X)\ p ] < oo 

(2 A>0,xe[0,A 1 / d ] d A,,VeC 

and for all A > 0 sup E[|£(x;"Pa U X)\ p ] < oo. 

iSK d ,AgC 

We implicitly assume for all Z > 0 that fj := £(x, X (T Bi(x)) has moments 
no larger than those of £. 

2.2. General central limit theorems. Under stabilization and moment 
conditions, we will show in Theorem 2.1 that E/r^ K and Var/1^ K have vol¬ 
ume order asymptotics and that the scaling limit of the finite-dimensional 
distributions of the renormalized measures A a mean zero Gaus¬ 
sian field. Theorem 2.1 is a special case of the upcoming Theorem 2.4 and 
applications of both are described in Section 3. 

By the convergence of finite-dimensional distributions of random signed 
measures fi n to those of a generalized Gaussian field we mean the conver¬ 
gence in distribution of the integrals / f dp n to the corresponding normal 
random variables for all test functions / € C(A). This is the usual func¬ 
tional analytic point of view where a measure is viewed as a continuous 
linear functional acting on continuous functions. Henceforth we say that 
measures converge to a Gaussian field if their finite-dimensional distribu¬ 
tions converge. 

For all r > 0, let 

V%):= E[£ 2 (0;7U)] 

+ / (E[£(0; V T Uy) ■ £(y; V T U 0)] — E[£(0; 'P r )]E[£(? / ; V' t )])t dy, 
JR d 

where V' T denotes an independent copy of V T . 

Theorem 2.1. (i) //£ is stabilizing and satisfies (2.2) for some p> 1, 

then for all f € C{A) 

(2.3) lim E ^ ;/i ^ = [ f( x )E[^(0;V K(x) )]K{x)dx, 

A— >oo A J A 

whereas if £ is polynomially stabilizing and satisfies (2.2) for p = 4, then 

(2.4) lim V ai [(/i/U<-c)] _ j f 2 (^v^{K{x))^{x)dx. 

A— >oo A Ja 

(ii) If £ is exponentially stabilizing and satisfies (2.2) for all p > 0, then 
A _1/2 /4 K converges as A —> oo to a Gaussian field with covariance kernel 
Ja fi{x)f 2 (x)Vt{K(x))K(x) dx. 
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Statistical applications often require the analog of Theorem 2.1 for mea¬ 
sures induced by exactly n i.i.d. points on A. This “de-Poissonized” version 
of Theorem 2.1 goes as follows. Let X^, i > 1, be i.i.d. with common den¬ 
sity k, X n := {Xi,..., X n }, and pi := Yn=\ CnPQ; X n )5 Xi the random “de- 
Poissonized” measures induced by k and £. To obtain the convergence of the 
finite-dimensional distributions of pi'-= pi~ IE pi we need some additional 
terminology and assumptions. 

Let X m ^ n be a point process consisting of m i.i.d. random variables nf! d X 
on n l / d A, where X has density k. For all X, let H(X) := H^(X) := 
and for all A > 0 let H\(X) := C\( x 'j %)■ For any finite X, let A x (X) : = 
H(X Dx) — H(X). Say that H satisfies the bounded moments condition for 
k (cf. [24]) if 

(2.5) sup sup sup E[A^(T m]n )] < oo. 

n x£rP/ d A m£[n/2,3n/2] 

If H satisfies the bounded moments condition for k then we will assume 
throughout that defined by Hj(X) := J2 x ex f{ x )H x 'i X), f £ C(A ), also 
satishes the bounded moments condition. This assumption is satisfied in all 
of our applications in Section 3. 

The next definition recalls a notion of stabilization for H introduced in 
[24], We are grateful to Mathew Penrose for pointing out that stabilization 
of H rather than that of £ is essential for the upcoming de-Poissonization 
methods of Section 6; this observation corrects an earlier version of our 
results. 

Definition 2.3. The functional H := H% is strongly stabilizing if for 
all r > 0, there exist a.s. finite random variables S (a radius of stabilization 
of H) and A^(r) such that with probability 1, 

(2.6) A^(r)=Ao((? T nB s (0))u4) 
for all finite 4cR d \ B${ 0). 

If H^ is strongly stabilizing, then we will assume throughout that f £ 
C(A), is also strongly stabilizing. This assumption is satisfied in all of our 
applications in Section 3. 

Let D^(t) := E[A^(r)] for all r > 0. The following de-Poissonized version 
of Theorem 2.1 shows that E[(/, pi)) and Var[(/, pi)], f £ C(A), have volume 
order fluctuations and that the scaling limit of the re-normalized measures 
n~ l / 2 pi is a mean zero Gaussian held. 

Theorem 2.2. (i) If f is stabilizing and satisfies (2.2) for some p> 1, 

then for all f £ C(Al), 

(2.7) lim MlM= f f( x )EmV K{x) )Mx)dx, 

n—>oo U J A 
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whereas if £ is polynomially stabilizing and satisfies (2.2) for p = 4, and if 
H is strongly stabilizing and satisfies the bounded moments condition for n, 
then for all f G C(A), 


( 2 . 8 ) 


lim 

n— kx> 


Var[(/,/4)] 


^ n 

= J f 2 (x)V^(n(x))n(x) dx — {^J f (x)( k(x))k(x) dx 


(ii) If £ is exponentially stabilizing and satisfies (2.2) /or allp> 0 and 
if H is strongly stabilizing, then n _ 1 ,/ 2 p| converges as n —> oo to a Gaussian 
field with covariance kernel 


(2.9) 


fi(x)f 2 (x)V^(K(x))K(x) dx 


— / fi(x)D^(n(x))K(x)dx / f 2 (x)D^(K(x))n(x)dx. 
Ja J a 

(iii) If the distribution of A^(k(X)) is nondegenerate, then 

( 2 . 10 ) 




n 


that is, the limiting Gaussian field is nondegenerate. 


Remarks, (i) Theorems 2.1 and 2.2 generalize existing central limit 
theorems in geometric probability (Heinrich and Molchanov [13], Malyshev 
[19], Penrose and Yukich [24] and Ivanoff [16]) in several ways: (a) they 
show asymptotic convergence of measures to a Gaussian held, thus also 
yielding asymptotic convergence of functionals to a limiting normal random 
variable, (b) they identify the limiting variance and covariance structure 
in terms of the underlying density of points, and (c) they do not assume 
spatial homogeneity of the underlying points. Theorem 2.2 implies that for 
fi, ■ ■ ■ ,fm € C(A), the random vector <</i, n _ 1 / 2 p|),..., {f m , n" 1 / 2 p|)) con¬ 
verges to a multivariate Gaussian random variable. 

(ii) Evaluating (2.4) and (2.8) is in general difficult. However, for some 
problems of geometric probability, for example, those involving functionals 
which count the number of pairs of points within a specified distance of 
one another, it is relatively simple to evaluate V^ and D% [ 6 ]. Moreover, a 
simplification of (2.4) and (2.8) occurs whenever £ is homogeneous of order 
7 , that is, whenever there is a constant 7 > 0 such that £ satisfies the relation 
£(ax; aX ) = a 7 £(x; X) for all positive scalars a and all finite point sets X 3 x. 
Homogeneity occurs naturally in many problems of geometric probability. 
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If £ is homogeneous of order 7 , then V^(t) = V^(l)r 27//rf , and D^{t) = 
D^( l)r — yielding 


( 2 . 11 ) 


lim 

A— xxd 


Var [(/^L)] 

A 



and 


( 2 . 12 ) 


lim 

n—>-c50 


Var[(/,p|)] 


n 


V^{1) jj 2 {x)K{x) { - d ~ 2 ^/ d dx 

- P ? (l)) 2 (/ A /(xM*)^" 7 ^ dx) 2 . 


If k is the uniform distribution on the unit cube, then by (2.12) 
lim Var ^([ Q,1 ] d) ] = V*(l) - (D«(l)) 2 , 

n—>00 77, 

which is strictly positive whenever A^(l) is nondegenerate. If £ is scale 
invariant, or homogeneous of order 0, then for any k with support A, 


i im AMAI = y£ ( i) _ ( 0 i (1) )7 

n—>oo ji 

showing that the limiting variance is not sensitive to the underlying density 
but depends only on the dimension. 

Still in the setting of general k, the inequality f A n(x)( d ~ 2 ' y ^ d dx > (f A k(x ^d-^/d dx)2 
implies that the right-hand side of (2.12) is strictly positive whenever A^(l) 
is nondegenerate. Moreover, (2.12) implies that when d = 2, k a density 
on A = [0,1] 2 , and 7 = 1, which would be the case for total edge length 
functionals of graphs on vertex sets in [ 0 , 1] 2 , then the limiting variance 
of equals 1^(1) — (D^(l)) 2 (f A k(x) 1 / 2 dx) 2 , which is minimized 

when the underlying density k is uniform. 

(iii) A comparison of (2.4) and (2.8) shows that Poissonization contributes 
extra randomness which shows up in the limiting variance (2.8). To show 
nondegeneracy of A^(k(X)), we need to appeal to the particular geometric 
structure of the underlying problem. This is done on a case by case basis and 
is already treated in many problems of interest [24]. The implicit finiteness 
of the right-hand side of (2.4) and (2.8) is made explicit in Section 4.3. 

(iv) Our method of proof actually yields (2.4) whenever (2.2) holds for 
some p > 2 and £ is exponentially stabilizing. This modification requires a 
small modification of Lemma 4.2. Also, if £ satisfies stabilization (Definition 
2.1) only when 0 < a < b, then Theorems 2.1 and 2.2 hold provided that k 
is bounded away from zero. 

(v) The condition m € [n/2,3n/2] in (2.5) is needed in order to achieve 
an efficient coupling between Poissonized and de-Poissonized measures. See 
Lemma 6.2 for details. 

(vi) Theorem 2.1 holds for arbitrary continuous K:M d —that is, k 
need not be a probability density. 
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2.3. Extensions of main results. 

2.3.1. Random measures induced by marked point processes. Theorems 
2.1 and 2.2 extend to random measures induced by marked point processes. 
Let (Abe a probability space of marks and let V TXV (resp. V KXU ) 
be a Poisson point process on x Af with intensity measure r x v (resp. 
k x u). We say that £ stabilizes if Definition 2.1 holds with Br(x) replaced 
by B^(x) := Br{x) x Ad, and X ranging over the finite subsets of (M d \ 
Br(x)) x M. Write f,(x;V T ) for f(x]V TXl/ ). 

Let Xi,i > 1, be i.i.d. marked random variables with common law dn x v. 
X- denotes the projection of Xi on M rf , X n := {Xi,...,X n }, and : = 
Ya=i £n(Xi m , X n )5x'. the associated marked random measures on M. d . Let 
A£(r) denote the marked version of (2.6), that is, 

/n -i o\ A $ (r):=A^(rxi/) 

1 ' j =A 0 ((P TXV ns^(o))ui) 

for all finite A C (M d \i? 5 ( 0 )) x AL Let V^{r) := V^(t x v) denote the marked 
version of V^(t). 

The analog of Theorem 2.2 for marked processes is as follows: 

Theorem 2.3. Let p^,n > 1, denote the marked measures defined above. 
Then (2.7)-(2.10) hold with D^(r) and V^(t) replaced by their respective 
marked versions E[A^(r x v)\ and V^(r x u). 

Remarks, (i) Theorem 2.3 generalizes Theorem 3.1 of [25], which es¬ 
tablishes a CLT for functionals of marked homogeneous samples. 

(ii) Applications of Theorem 2.3 to random sequential packing and spatial 
birth-growth models are discussed in Section 3.2. 

2.3.2. Random measures induced by nontranslation-invariant functionals. 
It takes just a little extra effort to use our general approach to prove CLTs 
for measures induced by nontranslation-invariant weights f. Although trans¬ 
lation invariance is often present in the measures (i)-(iii), we envision situ¬ 
ations where measures on do not enjoy translation invariance, as would 
be the case if the metric on M. d changes from point to point. 

Let £(y,x,X),y € be a family of measurable M-valued functions de¬ 
fined for all pairs (x,X), where X C M. d is finite and x is an element of A. In 
cases with x ^ X, we abbreviate the notation £(y;x, I’Ll {x}) to £(y;x, X). 
We assume for all y that £(y; x, X) is translation invariant in the pairs (x, X), 
that is, for all zGl 1 * and all pairs (x, X ), 

€{y;x,x) = £(y,x - z,x - z). 
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£(y;-,-) is a rule depending on y E M. d which assigns a real value to all 
pairs (x,X). We do not assume that £ is translation invariant in the triples 
(y,x,X). We define £(x;ff) := £(x\x,X) for all x E and for all A > 0, we 
set 

Z\(y,x,X):=£(y; X 1/d x,X 1/d X) and £a(x; X) := £(x; X 1/d x, \ 1/d X). 


We will consider limit theorems for the random measures 


/4 : = 5Z 

and 

n 

(2-14) 

2—1 


£ is said to be stabilizing if £(x;ff) = £(x;x,ft) stabilizes in the sense of 
Definition 2.1. If the rules £(y; •, •) are identical for all y E then £(x — 
z; ft — z) = £(#; ft) for all z E M d , and we reduce to the translation-invariant 
setting of Theorems 2.1 and 2.2. 

When £ is translation invariant, that is to say, when £(x;x,ft) =£(x,ft), 
stabilization guarantees that pair correlation functions for £ decay suitably 
fast enough with respect to the interpoint distance. It also guarantees that 
the pair correlation function with respect to nonhomogeneous samples can 
be closely approximated by the pair correlation function with respect to 
homogeneous samples. However, for nontranslation-invariant £, a suitable 
approximation of pair correlation functions is not possible without some 
continuity of £ with respect to its first argument. This motivates the follow¬ 
ing definition. 


Definition 2.4. The function £ is slowly varying in L q (abbreviated £ E 
SV(< 7 )) if for all t E (0, oc), any x E M cZ , and any compact set K containing 

0 : 

lim sup E[|£ a (x + x, V\ T ) - £,\(x;x,V\ T )\ g ] =0. 

A—>oc y£ K 

The following generalizes Theorem 2.1 to nontranslation-invariant £. For 
all x E and r > 0, let 

V*(x,r) :=E[£ 2 (x;x, 1P t )] 

+ [ (E[£(x;0,P t Ui/).£(s;i/,P t U0)] 
jR d 

-E[£(x;0,P T )]E[£(x;y,^)])rdj/. 
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Theorem 2.4. (i) //£ £ 5V(|) is stabilizing and satisfies (2.2) for some 

p > 1, then for all f £ C{A), 

(2.15) lim = I f( x )E[f(x]V K ( x) )]K,(x)dx, 

A— >oo A J A 

whereas iff is polynomially stabilizing and satisfies (2.2) forp = 4, then for 
all f £ C(A), 

(2.16) lim 1 = f f 2 (x)V^(x,K(x))K,(x)dx. 

A—>oo A J A 


(ii) If f £ SV(-) exponentially stabilizing and satisfies (2.2) /or all 
p > 0, f/ien A“ 1/2 /^ k converges as \ —+ oo to a Gaussian field with covariance 
kernel f A fi(x)f 2 (x)V^(x, k(x))k(x) dx. 

Letting H^-(X) = ffx £ xf(x-,X), we say that H is strongly stabilizing for 
all r > 0 and all x £ if there exist a.s. finite random variables S (a radius 
of stabilization of H) and A^(r,x) such that with probability 1, 

(2.17) A^(r,x) = A x ((P r nB s (x))uA) 

for all finite A C l d \ Bg (x). 

It is easy to check that the measures (2.14) satisfy the law of large num¬ 
bers in Theorem 2.2(i). The following de-Poissonized version of Theorem 2.4 
generalizes Theorem 2.2 and shows that the normalized versions of the mea¬ 
sures (2.14) converge to a Gaussian field as well. Put D^(x,r) := E[A^(x, t)]. 


Theorem 2.5. Let f £ SV(|). Assume that H is strongly stabilizing 
and satisfies the bounded moments condition for n. Then we have: 

(i) Iff is polynomially stabilizing and satisfies (2.2) forp = 4, then for 
all feC(A), 


(2.18) 


lim 

n—>oo 


Var[(/,4)] 


n 


= j f 2 (x)V^(x,n(x))K(x)dx— ^ J f(x)D^(x, k(x))k(x) dx 


(ii) Iff is exponentially stabilizing and satisfies (2.2) for allp> 0, then 
n ~ 1 ^ 2 Pn converges as n —> oo to a Gaussian field with covariance kernel 


(2.19) 


fl(x) f 2 (x)Vfix, k(x))k(x) dx 
- f fi(x)D^(x, k(x))k(x) dx f h(x)DHx,K(x)Mx)dx. 
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(iii) If the distribution of A^(X, k(X)) is nondegenerate, then 


lim 

n—xx) 


Var[/4(v4)] 

n 


> 0 . 


Remarks, (i) Formulas (2.18) and (2.19) are, in general, difficult to 
evaluate explicitly. However, in the context of statistics involving one-dimensional 
spacings, these formulas are readily evaluated [6], thus extending existing 
CLTs for sum functions of spacings. 

(ii) We have used the assumption £ £ SV { 4/3) only for technical con¬ 
venience and have not aimed to find the optimal choice of SV(q). Higher 
moment assumptions on £ will in general require £ £ SV ( q ) for smaller values 
of q (cf. Lemma 4.2). 


2.3.3. Random measures induced by graphs. Theorems 2.1 and 2.2 as¬ 
sume a special form when the random point measures are induced by graphs. 

We see this as follows. Let X be a locally finite point set and let G := G{X) 
be a graph on X. G is translation invariant if translation by y is a graph 
isomorphism from G(X) to G(y + X ) for all y £ M. d and all locally finite X. G 
is scale invariant if scalar multiplication by a induces a graph isomorphism 
from G{X) to G(aX) for all X and all a > 0. Given G and a vertex x £ X, 
let £(x]G(X)) be the set of edges incident to x (or for the Voronoi graph, 
the set of edges whose planar duals in the Delaunay graph are incident to 
x), and let |e| denote the length of an edge e. 

For any / £ B(a,b), let Vf iX denote Vf together with a point at x. 

Definition 2.5. G stabilizes if for all A € A' , 0 < a < b < oo, A > 0, and 
x £ A A, there exists an a.s. finite random variable R(x) := R(x, X,a,b, A) 

(a radius of stabilization ) defined on (Q,iF,P) such that for all / £ iF(a,b), 
with Supp / = AH, and all finite X C AH \ B R (x), we have 

£(x- G(Vf, x n B r (x) U X)) = £{x; G(V f , x n B R (x))). 

Given </>: R + —> R + , consider functionals of the type £^(x; X) := Jf e e£(x-,G(X)) ‘KM) 
such functionals could represent, for example, the total length of ^-weighted 
edges in G incident to x, the number of edges in G incident to x, or the 
number of edges in G less than some specified length. These functionals 
induce the point measures 

Y, ‘KMKz- 

xeX e£S(x;G(X)) 
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If G is polynomially stabilizing (resp. exponentially stabilizing), then so is 
for any fi. Given p > 1, say that £9 is L p bounded if 


( 2 . 20 ) 


sup sup sup E 
A>Oo; 6 a VdAXGC 


E 


eeS(x;G(V XK UA)) 


< oo. 


Let H^(X) denote the total mass of pG x - H^(X) is strongly stabilizing 
if for all r > 0 there exist a.s. finite random variables S := S{t) and A^'(r) 
such that with probability 1 

H$(V T u { 0 } n Bs( 0 ) \JA)~ H°(V T n b s ( o) u a) = a£(t) 

for all finite A C M. d \ Br( 0). 

Let X n := n l / d {X i,..., X n ), with X{ i.i.d. with density k. Write 


( 2 . 21 ) 


Pt,n : =J2 J2 <Ml e IMw- 

* =1 eGSi^/d-Xi^Xn)) 


The next result is the CLT counterpart to the main result of [26] and 
follows immediately from Theorem 2.2. There is obviously a Poisson version 
of Theorem 2.6, but we do not state it here. Put D^(t) := E[A(p(r)] and let 
i/jp be the function W when £ := is defined as above. 


Theorem 2.6. Assume that the graph G is translation and scale in¬ 
variant. Let Xi,i > 1, be i.i.d. with density n. Assume that is strogly 
stabilizing and satisfies the bounded moments condition. 

(i) If satisfies (2.20) forp = 4, if G is polynomially stabilizing, then 
for all t > 0, 


lim 

A—>oo 


Var[frG(Ai/ d ( P AT n[0,l] d ))] TrG 


= V* (r) ■ T 


and for all f E C{A), 

Var [(/,P?,n)] 


lim 

( 2 . 22 ) rwo ° 


n 


= J a f 2 { x ) v ^(k(x))k( x ) doc - (^jj{x)D^{K.{x))K{x)dx S j . 

(ii) If G is exponentially stabilizing, if satisfies (2.20) for all p > 0, 

then as n —► oo, n _1 / 2 p^ n converges to a Gaussian field with covariance 
kernel 


fl(x)f 2 (x)V$(K(x))K(x) dx 


- I fi{x)D$(n(x))K(x)dx I f 2 (x)D < £(K(x))K(x)dx. 


G( 


A 


(2.23) 
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Remarks, (i) If 4>(x) = x p , p > 0, then the integrals in (2.22) and (2.23) 
can be simplified using the identities V/ g (t) = V9(l )r~ 2p / d and D^(r) = 
D^\)r-P/ d . 

(ii) We may generalize Theorem 2.6 to treat nontranslation-invariant £. 
For example, let ^{x\x,X) be a functional which assigns to a point x in 
the graph G(X) a value which depends on the point iGK' 1 (e.g., the value 
may depend upon the local metric structure at x). Such functionals are 
not translation invariant in the triples (x',x,X). By applying an appropriate 
uniformization to curved surfaces, we can fit functionals on such surfaces 
into our set-up of nontranslation-invariant functionals of point processes on 
R d . This yields, for example, CLTs for functionals of graphs defined over 
curved surfaces, in particular functionals of Voronoi diagrams over surfaces 

[ 17 ]. 

3. Applications. Theorems 2.1-2.6 can be applied to point measures in¬ 
duced by random graphs, packing processes and germ-grain models. This 
extends previous results [2, 21, 24, 26], [5, 7, 8, 9, 11, 15, 25, 27] and 
[12, 13, 22, 26] to the weak limit setting as well as to the setting of in¬ 
teraction processes over nonhomogeneous point fields. We do not provide an 
encyclopedic treatment of applications and anticipate applications to other 
interaction processes on R d , including measures induced by continuum per¬ 
colation models. The methods described here can be modified to extend and 
generalize the central limit theory for classical spacings and (^-divergences; 
in this setting the functions V * and may be determined explicitly, al¬ 
lowing us to compute the limiting variance explicitly as a function of the 
underlying density of points. We refer to [6] for complete details. 

Throughout, we will often show the exponentially stabilizing condition by 
appealing to results of [24, 26], which involves a slightly stronger definition 
of stabilization. 

3.1. Random graphs. We limit discussion to random graphs on R d with 
the usual Euclidean metric, but since translation invariance of £ is not as¬ 
sumed, many results hold if the graphs are defined on curved spaces. Our 
discussion parallels that in [26]. We say that <f> has polynomial growth if there 
exists a < oo such that (j)(x) <(7(1 + x a ) for all x G R + . 

3.1.1. k-nearest neighbors graphs. Let A: be a positive integer. Given a 
locally finite point set X C R rf , the A;-nearest neighbors (undirected) graph 
on X, denoted NG(T), is the graph with vertex set X obtained by includ¬ 
ing {x,y} as an edge whenever y is one of the k nearest neighbors of x 
and/or x is one of the k nearest neighbors of y. The /c-nearest neighbors 
(directed) graph on X, denoted NG^T), is the graph with vertex set X 
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obtained by placing a directed edge between each point and its k nearest 
neighbors, /c-nearest neighbors graphs are translation and scale invariant. 
Given a binomial sample X\,... ,X n of i.i.d. random variables with density 
ft, define the induced point measures pj^ G and p]^ G as in (2.21). 

Theorem 3.1. The random measures p^ G and p^f G , n > 1, satisfy 
(2.22), (2.23) if 4> has polynomial growth, Supp ft G A, and k is bounded 
away from infinity and zero on its support. 


If we set (f> G ( |e|) = |e|/2, then we obtain a CLT for the total edge length of 
the fc-nearest neighbors graph on the nonhomogeneous point set X n when¬ 
ever Supp ft £ A and ft is bounded away from infinity and zero. This general¬ 
izes existing CLTs [2, 24] which only show CLTs for nearest neighbor graphs 
on homogeneous point sets. The convergence to a Gaussian limit (2.23) is 
new. 

Still more generally, if <^ G (|e|) = |e| p /2, p > 0, then Theorem 3.1 yields a 
CLT for the pth. power-weighted total edge length of the fc-nearest neighbors 
graph on X n when Supp ft £ A and ft is bounded away from infinity and 
zero. That is, there are constants P NG (1) and Z1 NG (1) such that 


Var [(/,P%n)] 


lim 

n—xx) 


(3.1) 


= P NG (1) / / 2 (x)ft(x) (d “ 27)/d dx 
JA 

- (H NG (1)) 2 f(x)K(x)^/ d dx 


and n 1//2 p^ G converges as n - 


oo to a Gaussian field with covariance kernel 


(3.2) 


P NG (1) / fi(x)f 2 (x)K(x)^- 2 ^/ d dx 
JA 

~(L> ng ( 1)) 2 / /i(x)ft(x) (d - 7)/d dx / f 2 (x)K(x) {d ~^ /d dx. 
JA JA 


Another application of Theorem 3.1 goes as follows. Fix t > 0. Let <p G (\e\) 
be either 0 or 1 depending on whether the length |e| of the edge e is bounded 
by t or not. Then (2.23) gives a CLT for the empirical distribution function 
of the rescaled lengths of the edges in the /c-nearest neighbors graph on X n . 


Proof of Theorem 3.1. The proof is straightforward and essentially 
follows from existing arguments in [24] and [26] . For completeness we sketch 
the proof when G{X) denotes NG(ff); similar arguments apply when G(X) 
denotes NG^T). It will suffice to apply Theorem 2.6 and to show that 
NG stabilizes on elements of J-(a,b ) when a > 0 [recall Remark (iv) after 
Theorem 2.2]. Let / G J-(a, b ) be arbitrary, where 0 < a < b < oo. As shown in 
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Lemma 6.1 of [24] (even though the definition of stabilization there is slightly 
different), the set of edges incident to x in NG(7 y f !X ) is unaffected by the 
addition or removal of points outside a ball of random almost surely finite 
radius 4 R, that is, the graph G(X) = NG(T) is stabilizing. Moreover, R is 
constructed as follows [24]. For each t > 0 construct six disjoint equilateral 
triangles Tj(t), 1 < j < 6, such that x is a vertex of each triangle and such 
that each triangle has edge length t. Then R is the minimum t such that each 
triangle Tj(t), 1 < j < 6, contains at least k + 1 points from Vf_ x . Since / is 
bounded away from zero, elementary properties of the Poisson point process 
give the desired exponential decay of R and thus 4 R decays exponentially 
as well. We verify the moments condition (2.20) as in the proof of Theorem 
2.4 of [26]. Strong stabilization of H is given by Lemma 6.1 of [24] and the 
bounded moments condition for H is as in Lemma 6.2 of [24]. The positivity 
of the limiting variance is given by Lemma 6.3 of [24]. □ 

3.1.2. Voronoi and Delaunay graphs. Given a locally finite set X C U d 
and x € X, the locus of points closer to x than to any other point in X is 
called the Voronoi cell centered at x. The graph on vertex set X in which 
each pair of adjacent cell centers is connected by an edge is called the De¬ 
launay graph on ff; if d = 2, then the planar dual graph consisting of all 
boundaries of Voronoi cells is called the Voronoi graph generated by X. 
Edges of the Voronoi graph can be finite or infinite. Let DEL(T) [resp. 
VOR(ff)] denote the collection of edges in the Delaunay graph (resp. the 
Voronoi graph) on X. The Voronoi and Delaunay graphs are clearly scale 
and translation invariant. Define the induced point measures p^® R and p RRL 
as in (2.21). 

Theorem 3.2. Let d = 2. The random measures /5^® R and p RRL satisfy 
(2.22), (2.23) if has polynomial growth with 4>(oo) = 0, Supp k € A and if 
k is bounded away from infinity and zero on its support. 

The limits (2.22), (2.23) extend the results of Penrose and Yukich [24] and 
Avram and Bertsimas [2], who consider CLTs for the total edge length of 
Voronoi graphs over homogeneous samples. The convergence in distribution 
to a Gaussian limit (2.23) is new. Clearly, the analogs of (3.1) and (3.2) hold 
for the measures p^® R and p RRL . 

Proof of Theorem 3.2. We will apply Theorem 2.6. The moments 
condition (2.20) is shown in Theorem 2.5 of [26]. We can verify as in [26] 
that G(X) is stabilizing. Let / £ T{a,b) be arbitrary, with 0 < a < b < oo. 
We will show that the Voronoi cell centered at x with respect to Vf x is 
unaffected by changes beyond a random but a.s. finite distance R. from x. 
We only need to show R has exponentially decreasing tails. This is done 
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in a manner similar to that for the fc-nearest neighbors graph. For each 
t > 0 construct 12 disjoint isosceles triangles Tj(t), 1 < j < 12, such that x 
is a vertex of each triangle, such that each triangle has two edges of length 
t, where Tj(t) C Tj(u ) whenever t < u and where U* > o U}=i = K 2 - If 
R is the minimum t such that each triangle Tj(t), 1 < j < 12, contains at 
least one point from Vf tX , then 3 R is a radius of stabilization (page 1037 
of [24]). Since / is bounded away from zero, elementary properties of the 
Poisson point process give the desired exponential decay of 3 R. We can 
verify the bounded moments condition on H as in Lemma 8.1 of [24]. Strong 
stabilization of H is proved in Section 8 of [24] . The positivity of the limiting 
variance is given by Lemma 8.2 of [24], □ 

3.1.3. Sphere of influence graph. Given a locally finite set X C M d , the 
sphere of influence graph SIG(Af) is a graph with vertex set X, constructed 
as follows: for each x € X let B(x) be a ball around x with radius equal to 
m i n y6A\{z}{|y — a:|}- Then B{x) is called the sphere of influence of x. Draw 
an edge between x and y iff the balls B(x) and B{y ) overlap. The collection 
of such edges is the sphere of influence graph (SIG) on X and is denoted 
by SIG(T’). It is clearly translation and scale invariant. Define the induced 
point measure p ^ as in (2.21). 

In Section 7 of [24], CLTs are proved for the total edge length, the number 
of components, and the number of vertices of fixed degree of SIG when the 
underlying sample is uniform. The following extends these results to nonuni¬ 
form samples and also shows weak convergence of the associated measures. 
We also obtain a CLT and variance asymptotics for the total number of 
edges in the SIG on nonuniform samples, extending results of [14]. 

Theorem 3.3. The random measures p® 1 ® satisfy (2.22), (2.23) if f 
has polynomial growth, SuppK £ A and k is bounded away from infinity and 
zero on its support. 

Proof. We will apply Theorem 2.6 again. Let / b) be arbitrary, 

0 < a < 6 < oo. As shown in [26], G(Vf) has moments of all orders and is 
stabilizing, so we only need to show exponential stabilization. However, this 
follows from the analysis of SIG in [24], Consider an infinite cone C with 
its vertex at x, subtending an angle of 7r/6 radians. As in [24], let T be the 
distance from x to its closest neighbor in Vf )X D C , and if Y is the point in 
C I^B§t{x) closest to x, then note (page 1030 of [24]) that the configuration 
of points outside B 3 \ Y \(x) has no effect on the set of points in C connected to 
x. Thus, the radius of stabilization R equals the maximum of m i.i.d. copies 
of 3|T|, where m is the minimum number of cones C),..., C m congruent 
to C , each with vertex at x, whose union is M. d . It is easy to check that R 
has exponential tails. Thus G{X ) is exponentially stabilizing. The bounded 


18 


YU. BARYSHNIKOV AND J. E. YUKICH 


moments condition on H is as in Lemma 7.2 of [24] and strong stabilization 
of H is as in Lemma 7.1 of [24], The positivity of the limiting variance is 
given by Theorem 7.2 of [24]. □ 

3.2. Random packing. We will use Theorem 2.3 to extend earlier results 
on random sequential packing [5, 7, 8, 9, 11, 15, 25, 27] to cases of nonhomo- 
geneous input as well as to show the weak convergence of packing measures 
induced by Poisson and fixed input. 

3.2.1. RSA packing. The following prototypical random sequential pack¬ 
ing model is of considerable scientific interest. 

Let B ni \,B n ^ 2 , ■ ■ ■, B n ^ n be a sequence of d-dimensional balls of volume n -1 
whose centers are i.i.d. random d- vectors X \,..., X n with probability density 
function k : A —> [0, oo). Without loss of generality, assume that the balls are 
sequenced in the order determined by marks (time coordinates) in [0,1]. 
Let the first ball B n ^i be packed, and recursively for i = 2,3,..., N, let the 
Ah ball B Ui i be packed iff B n> i does not overlap any ball in B n> i,... ,B Ht i-\ 
which has already been packed. If not packed, the ith ball is discarded. The 
collection of centers of accepted balls induces a point measure on A, denoted 
H n -i. We call this the random sequential packing measure induced by balls 
(of volume n -1 ) with centers arising from k. 

Packing models of this type arise in diverse disciplines, including physi¬ 
cal, chemical and biological processes. In statistical mechanics, this model 
describes the irreversible deposition of colloidal particles or proteins onto a 
substrate. In this context, the model described above is known as the RSA 
model for hard spheres on a continuum substrate. When the ball centers 
belong to a stationary Poisson point process on W l , then this model is the 
Matern hard-core process (page 163 of [29]). When the ball centers belong 
to a bounded region of M rf , then this model is known in spatial statistics as 
the simple sequential inhibition model (page 308 of [30]). 

The vast scientific literature on versions of RSA models (see [25] for ref¬ 
erences) contains an abundance of experimental results, but few rigorous 
mathematical results. In d= 1, Renyi [27] and Dvoretzky and Robbins [11] 
established LLNs and CLTs, respectively, for the total number of accepted 
balls. Coffman, Flatto, Jelenkovic and Poonen ([9], equation (2), Theorems 
13 and 14) determine explicit formulae for some of the limiting constants 
in the LLN and CLT, but restrict attention to d = 1. In d > 1, Penrose and 
Yukich [25] establish the asymptotic normality of the number of accepted 
balls when the spatial distribution is uniform and also show [26] an LLN for 
the number of accepted balls when the spatial distribution is nonuniform. 
Baryshnikov and Yukich [5] establish weak convergence of the sequential 
packing measures in d> 1, but only for homogeneous Poisson input. Here 
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we will use our general result for marked processes, Theorem 2.3, to estab¬ 
lish convergence of the variance and also weak convergence in the case of 
nonhomogeneous input in d > 1. We will follow the set-up of [26] . 

For any finite point set X C A, assume the points x € X have time coor¬ 
dinates which are independent and uniformly distributed over the interval 
[0,1]. Assume balls of volume n” 1 are centered at the points of X and 
arrive sequentially in an order determined by the time coordinates, and as¬ 
sume as before that each ball is packed or discarded according to whether 
or not it overlaps a previously packed ball. Let £(x\X) be either 1 or 0 
depending on whether the ball centered at x is packed or discarded. Let 
£ n (x;X) = £(n 1//rf x; nf! d X\ where n}^ d x denotes scalar multiplication of x 
and not the mark associated with x and where balls centered at points of 
n l / d X have volume 1. Let H(X) -=Y1xgx £(x\X) be the total number of 
balls packed. The random measure 

n 

&-=Y,Zn{Xi\{Xi}U)8x i 

i =1 

coincides with n n -i. 

Straightforward modifications of [25, 4] show that £ is exponentially sta¬ 
bilizing. The strict positivity of V^(t) is shown in Theorem 1.2 of [25]. Since 
£ is bounded it satisfies the moments condition (2.2). By Section 5 of [25], H 
satisfies the bounded moments condition and strong stabilization. Therefore, 
Theorem 2.3 yields the following CLT. 

Theorem 3.4. The random measures //^,n > 1, satisfy (2.8) and (2.9). 

Theorem 3.4 shows asymptotic normality of the total number of accepted 
balls and generalizes [4, 25] to the case of nonhomogeneous input. 

3.2.2. Spatial birth-growth models. Consider the following spatial birth- 
growth model in M rf . Seeds are born at random locations Xi £ M. d at times 
T t , i = 1,2,..., according to a unit intensity homogeneous spatial temporal 
Poisson point process 'L := {(W,Tj) £ M. d x [0, oo)}. When a seed is born, 
it forms a cell by growing radially in all directions with a constant speed 
v > 0. Whenever one growing cell touches another, it stops growing in that 
direction. Initially the seed takes the form of a ball of radius pi > 0 centered 
at Xi. If a seed appears at Xi and if the ball centered at Xi with radius pi 
overlaps any of the existing cells, then the seed is discarded. 

In the special case when the growth rate v = 0 and pi is constant, this 
model reduces to the RSA packing model. In the alternative special case 
where all initial radii are zero a.s., the model is known as the Johnson- 
Mehl model, originally studied in model crystal growth, and is described in 
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Stoyan, Kendall and Mecke [29]. Chiu and Quine [8] show that the number 
of seeds accepted inside a cube Q\ of volume A by time t satisfies a CLT, 
but apart from numerical considerations, their arguments do not preclude 
the possibility of limiting normal random variable with zero variance [7]. 
Penrose and Yukich [25] consider a modification of this model in which all 
seeds outside Q\ are automatically rejected, while the rules for seeds inside 
Qx x [0, oo) are as above. They establish a CLT for this model and show 
that the limiting variance is strictly positive (page 295 of [25]), thus implying 
that the CLT of [8] is nondegenerate. 

If seeds are born at random locations AQ £ A, it is natural to study the 
spatial distribution of accepted seeds. As far as we know, this problem has 
not been investigated. We may use Theorem 2.3 to establish the weak con¬ 
vergence of the random measure induced by the locations of the accepted 
seeds. 

For any finite point set X C A, assume the points x £ X have i.i.d. marks 
over [0,1]. A mark at x £ X represents the arrival time of a seed at x. 
Assume that the seeds are centered at the points of X and that they arrive 
sequentially in an order determined by the associated marks, and assume 
that each seed is accepted or rejected according to the rules above. Let 
£(x;X) be either 1 or 0 according to whether the seed centered at x is 
accepted or not. H ( X ) := Y^xex £0*3 X) is the total number of seeds accepted 
and A^(r) is as in (2.7). 

As with RSA packing, let X \,..., X n be i.i.d. random variables with den¬ 
sity k on A and with marks in [0,1]. The random measure 

n 

0 $ :=£&(*; {*<}£=i)(S* 

i— 1 

is the scaled spatial birth-growth measure on A induced by Xi ,..., X n . The 
next result, a consequence of Theorem 2.3, shows that the spatial birth- 
growth measures converge to a Gaussian field. 

Theorem 3.5. The random measures <r|,n> 1, satisfy (2.8) and (2.9). 

Theorem 3.5 generalizes [8] and extends [25] to the case of nonhomoge- 
neous input. 

3.2.3. Related packing models, (a) Theorem 3.4 extends to more general 
versions of the prototypical packing model. By following the stabilization 
analysis of [25], one can develop asymptotics in the finite input setting for 
the number of packed balls in the following general models: (i) models with 
balls replaced by particles of random size/shape/charge, (ii) cooperative se¬ 
quential adsorption models and (iii) ballistic deposition models (see [25] for 
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a complete description of these models). In each case, Theorem 2.3 yields 
weak convergence to a Gaussian limit of the random packing measures asso¬ 
ciated with the centers of the packed balls, whenever the balls have a density 
k : A —* [0, oo). 

(b) The above packing models describe convergence of measures arising 
as a result of dependently thinning a Poisson point process. Related ways of 
thinning processes include the annihilating process, described as follows. A 
clock is attached to each point (particle) in the process; when the clock for 
a chosen particle rings, then if the particle has itself not been annihilated, it 
annihilates its neighbors within a fixed radius. Clearly, once a particle is free 
from occupied neighboring sites, it remains there undisturbed and is fixed 
for all time. Thus in any finite region the process is unchanging after a finite 
time. This models the thinning of seedlings [31] and the resulting random 
point measure satisfies the CLT in Theorem 3.4. 

3.3. Germ-grain models. Germ-grain models form a central part of stochas¬ 
tic geometry and spatial statistics [12, 22]. Here we consider the limit theory 
of functionals and measures associated with germ-grain models. Such mod¬ 
els fall within the scope of the general set-up of Heinrich and Molchanov [13], 
who were the first to develop a general limit theory for random measures 
induced by translation-invariant germ-grain models. 

Let T t ,i > 1, be i.i.d. bounded random variables defined on (H,JF, P), 
independent of the i.i.d. random variables AQ, i > 1, which are also defined on 
and which have density k. For simplicity, consider random grains 
having the representation AQ + n~ l / d BT, (AT,) and consider the random set 

n 

~ n :=U(AQ + n- 1/d i? Ti (0)). 

2—1 

When the AQ,i > 1, are the realization of a Poisson point process, the set 
S n is a scale-changed Boolean model in the sense of Hall [12], pages 141 and 
233. Heinrich and Molchanov [13] exploit the translation invariance of such a 
model and establish a central limit theorem for the associated measures. For 
translation-invariant models, Heinrich and Molchanov [13] establish CLTs 
without assuming boundedness of T % . 

For all u € M d , let T(u ) be a random variable with a distribution equal 
to that of T\. For all x G M d and all point sets X C R rf , let V(x,X) be the 
Voronoi cell around x with respect to X. Given x G M d , let L(x,X ) denote 
the Lebesgue measure of the intersection of the random set U uex ^T(u)( u ) 
and V ( x , X). 

The volume measure induced by E n is 

n 

L n(Xi] X n )6 Xi 

i— 1 
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and the total volume of n 1 ^ d r. n is given by H L (n l / d X n ) := -^npQ; <T n ). 

Since T is bounded it follows that L is exponentially stabilizing and 
that H l is strongly stabilizing. Moreover, since the functional L is bounded 
by the volume of a Voronoi cell, it is clear from Section 3.1.2 that L satisfies 
the moment condition (2.2) for all p > 0 and that H L satisfies the bounded 
moments condition (2.5) for any k € 0 < a < b < oo. 

Therefore, for germ-grain models r. n given above we have thus proved: 


Theorem 3.6. Let the density k be bounded away from, infinity and 
zero. 

(i) For all f € C(A) 

lim Var [if, fin)] 

(3-3) — n 

= J f 2 (x)V L (x, k,(x))k(x) dx — ( J f(x)D L (x, k(x))k(x) dx) . 

(ii) As n —> oo, converges to a Gaussian field with covariance 

kernel 


(3.4) 


fl(x)f 2 {x)V L (x , k(x))k(x) dx 


- I h{x)D L (x,K{x))K{x)dx [ f 2 (x)D L (x,K(x))K(x)dx. 
Ja j a 


Remarks, (i) We have confined attention to one of the simplest germ- 
grain models. Instead of balls Bt, one could assume that the grains have 
some distribution on the space of convex subsets of M d . We have also limited 
our discussion to volume functionals, but it should be clear that the approach 
above readily extends to other spatial statistics, including total curvature. 

(ii) Theorem 3.6 shows that volume functionals satisfy a CLT over nonuni¬ 
form point sets, adding to results of [13] and [12], Chapter 3.4, involving the 
vacancy functional for germ-grain models. 

(iii) The LLN counterpart of Theorem 3.6 is established in [26] and is not 
reproduced here. 


4. Proof of variance convergence (Poisson case). 

4.1. Correlation functions. The proof of (2.16) uses the objective method [1] 
together with correlation functions. To illustrate the method, it is instructive 
to first prove the limit (2.3). Recall that for all / G C(A), A” 1 E[(/,/i| K )] = 
f A f(x)K[f\(x-, V\ k )]k(x) dx. The key observation lying at the heart of the 
objective method is that for any point x € A distant at least (K/ A) 1//rf from 
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dA, K large, £,\(x: V\ K ) is well approximated by the candidate limiting ran¬ 
dom variable £(a:; V K r x )) in the sense that as A —» oo 

( 4 . 1 ) < \Mx{x;V X k ) ~^\(x;V X k(x))\ 

+ Ie^aOg^a+x)) -E£(x;P k{3 .))| - > o, 

where V XK t x \ is a Poisson point process on W l with intensity Xk(x) coupled 
to V\ K as in the upcoming coupling (4.2). 

Indeed, to prove (4.1), for any point x £ A distant at least (K/X) l / d from 
dA , consider the event Fk,x(x) that the radius of stabilization R(X l ^ d x) 
at \ 1 / d x with respect to V K ^ is less than K and that V\ K = 'Pak(x) on 
b {k/\)V d{x). Then 

P[F K ,x(x) c } < P[R(X l / d x) > K\ + A f \k{ V ) - «(®)| dy. 

By stabilization, we have P[R(X l t d x) > K\ < e uniformly in A by choosing 
K large enough. For such K , the Lebesgue point property of x shows that 
the second term above can be made arbitrarily small for large A and thus 
P[ f k,\( x )\ < 2e for lar g e Bound |E£ A (x; V Xk ) - M\(x',P X k(x))\ by 

|e[(£a(*;Pa«) -Ca(®;^a«(x))) ■ 1 f ka (x)}\ 

+ |e[(6(^;^a k ) - ^x(x\V X k { X ))) ■ 1 /•k,a(.t)'-]I- 

The first term vanishes by the definition of F^ x (x). The second term is 
bounded by a multiple of e by combining Holder’s inequality, the assumed 
1 + 5 moment condition and P[F^ x (x)\ < 2s. Thus, for any point x € A 
distant at least ( K/\) l / d from dA , the first term on the right-hand side of 

(4.1) goes to zero and the second term goes to zero by stabilization. 

The proof of the variance convergence (2.16) is more involved and requires 
some extra terminology. Let V Xk be a Poisson point process equidistributed 
with and independent of V Xk , that is, V' Xn is a copy of V\ K . For all A € M + 
and x € A we introduce two auxiliary homogeneous (independent) Poisson 
point processes V Xr j x ) and 7^ / ) defined on (Q,F,P) such that: 

(i) V\ K ( X ) and P Xk ^ have constant intensity on A equal to Xk(x), 

(ii) V\ K and V Xk ( x ) are coupled in the sense that for any Borel subset 
B c A, 

(4.2) P[V Xk (B)^V Mx) (B)]<X I \K(y)- K (x)\dy, 

J B 

and the same is true for V Xn and P Xk ( x ) ■ 
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The proof of the variance convergence (2.16) approximates the correla¬ 
tions of £,\{x\V\ «), x £ M rf , by those of ix(x',V xk,( x )) ■> x £ M d . Thus, for all 
x £R d define 

q x (x) :=E[$,l(x;Vx K )] and q x (x) := E[^(s; Vx K {x))], 
as well as the pair correlation function 

cx(x , y ) := E[£ A (x; x, Vxk U y)£\(y; V , Vxk U x)] 

- E[£a(x; X, 7 ? a k )]E[ 6(y; y, V'x K )] , x, y £ R d , x + y. 

Abbreviating notation throughout and writing £a(x;"Pak U y)£\(y; 'Pxk U 
x) for ^a(x; x, Vxk U y)£x(y, y, 'Pxk u x), we also have the pair correlation 
functions in the homogeneous intensity case: 

ca(x, y) := E[£a(x; V Xk ( x ) U y)£\(y; V\ K {x) u x) 

- Z\(x-,Vx K ( x ))£\(y,Vx K(x) )\, x^y, 

and 

^(x,y) :=eK*(®;^a«(*) uy)^(y;^«w ux) 
-?A^;^A K (a ; ))?A(y;^L(a : ))]> 

Here we employ the notation £ x (z',X) for ^a(x; 2 ,A). Clearly, the corre¬ 
lations cx{x,y) and c^(x,y) are not symmetric in x and y, unlike cx- When 
A = 1 we write simply q(x) and c(x,y) for q\{x) and ci(x,y), respectively, 
and similarly for q, c and c x . Denote the integral of f £ C{A ) with respect 
to a Borel measure y on R d by (f,y)- Now 

(4.3) A -1 Var[(/, y{ K )] = A(/ 0 f, E[y{ K 0 y{ K - /x| K 0 /i'jfj), 

where / 0 / denotes the function /(x)/(y) on the product A x A, y®v 

stands for the product measure on A x A and y'\ K is just an independent 

copy of y\ K . Considering the diagonal and off-diagonal terms, we may rewrite 
the integral (4.3) in terms of correlation functions 

A _1 Var [(/,i4k>] 

= A [[ f(y)f(x)cx(x,y)n(x)K(y)dxdy+ [ f(x) 2 q x (x)^x)dx 
JJaxA Ja 


= f k(x)/(x) 
Ja 


f( x )qx(x) + A / f(y)cx{x,y)n{y)dy 
Ja 


dx. 


To show the desired asymptotics (2.16), we will first show for all x £ A 
not too close to the boundary of A, that as A —> oo 


(4.4) 


f{x)qx{x) + A / f{y)cx(x,y)n(y)dy 


f(x)q(x)+ f (x)c x (x, y)n(x) dy 
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Note that the bracketed expression in (4.4) is f(x)V^(x, k(x)). 

4.2. Properties of correlation functions. Showing the limit (4.4) requires 
some properties of correlation functions. Using the definitions and the trans¬ 
lation invariance of £(y; x, X) in the pairs (x, X), it is easy to verify that for 
all x,y € A: 

q x (x) = q(x) and c x (x, x + y) = c x (x, x + X l ^ d y). 

Also, if (2.2) holds for p = 2, then we have the following uniform bounds: 

(4.5) sup [max[q x (x),qx{x),cx(x,y),c x (x,y),c x x (x,y)]]<oo. 

x,y£ R d ,A>0 

Our next fact provides some crucial decay properties. Here and elsewhere 
C denotes a constant whose value may change from line to line. 

Lemma 4.1. Under the moment condition (2.2) with p = 4, we have 
{m&x(\cx{x,y)\,\cx(x,y)\,\c x x(x,y)\)\ < C{r(X 1/d \x - y\/2)) 1/2 . 

Proof. We prove only 

\c x (x,y)\<C(r(X 1/d \x -y\/2)) 1/2 

since the proof of the other two inequalities is identical. Let R x := R(X 1 ^ d x) 
and R y := R(X 1 ^ d y) be the radii of stabilization of £ for X l / d x and A l ^ d y, 
respectively. Let E := E XiVt x denote the event that R( A 1//rf x) and R(X 1 ^ d y) 
are both less than 5 := X 1 / d \x — y\/2 and note that P[E C ] < Cr(X 1 ^ d \x — y\/2). 
On E the stabilization balls Br x (A 1 ^®) and Br v (A l ^ d y) do not intersect and 
thus on E we have £ A (x; V Xk U y) = £ A (x; V Xt f) and £ A (y; V Xk Ux) = £ A (y; V Xk ) 
showing that 

|E[£ a (t; Vxk u y)£\(y; V Xk u x)] 

- E[£ a (x; V Xk 0 B s (\ 1/d x))£ x (y, V Xk 0 B s (X 1/d y))}\ < CP[E C } 

by Holder’s inequality and the moment condition (2.2) with p = 4. Using 
independence in the second expectation and the bound 

mU^Vx^BsiX^x))] -E[£ a (s;P A k )]| < CP[E C }. 

we are done. □ 

The next lemma shows that q x and c A are closely approximated by their 
“uniform versions” q x and c^, respectively. Compactness of A and the con¬ 
tinuity of k and / imply uniform continuity, so we fix moduli of continuity 
t K ,tf :M + -+M + such that for any x,y G A : \x — y\ <5, \k(x) — n(y)\ < t K (5 ) 
and \f(x)-f(y)\ <t f (5). 
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Lemma 4.2. Assume that £ £ SV(^) satisfies the moment condition 
(2.2) for p = 4 and is polynomially stabilizing. Then there exists a func¬ 
tion e: —► M + , decreasing to 0, and a function 5: M + —> M + , increasing to 
oo, such that 5 /A —> 0 and 

(i) Vx E A distant at least (<5/A) 1//rf from dA, \q x (x) — q X (x)\ < e(A), 

(ii) V x,y E A, each distant at least (<5/A) 1//rf from dA, 

5(\)\c x (x,y) - c x x (x,y)\ < e(A), 

(iii) as A —> oo, the function 6 satisfies 

S(X)t f ((6(X)/X) 1/d ) - 0; «5(A)f K ((5(A)/A) 1/d ) - 0. 

Proof. It is clear that one can always find a function 6 := 5(A) —> oo 
as A —> oo such that 6/X —► 0 and implication (iii) holds, and even more 
strongly, that 

(4.6) S(X) 2 t K ((S(X)/X) 1 ^ d )—*0 as A —> oo. 

Fix such 5. 

For any A > 0 and x E A, x distant at least (5/A) 1,/d from dA, consider 
the event Q x := Q x ,\,s that the radius of stabilization R(X l / d x) of £ is less 
than 6, and that V Xk = V Xk ( x ) on B^ /X)1/d (x), that is, X l/d V\ K = X 1/d V XK (x) 
on B 5 \/d{X 1 ^ d x). By polynomial stabilization, by definition of t K , as well as 
by the coupling estimate (4.2), the probability of the complement of is 

(4.7) P[n c x ] < L0 d 5t K ((c 5/A) 1/d ) + oi<5-“ 2 , 

where a\ > 0 and 02 > 2 are constants and ca d is the volume of the unit ball 
in M. d . Thus P[Xl x ] tends to 0 as A —> oo by (4.6). To show (i) we need to 
show \E[fl(x;V XK ) - f x (x-,V x ^ x) )}\ -> 0 as A->oo: 

\ne x (x-,v XK )-e x (x-,v XK{x) )}\ 

< e [|£a (x-,V X k ) - fl(x;V XK(x) )\ln x \ 

The first term vanishes by the definition of Q x and the definition of £ A (x; V Xk / x \) 
£(A 1//d x; P K (x)). Holder’s inequality, the moment condition (2.2) with p = 4, 
and the bound (4.7) show that the second term vanishes as A —> oo. This 
proves (i). 

For the proof of (ii) it suffices to show that there exists a function 5(A) —► 
oo such that both 


(4.8) 


<5(A)|c A (x,y) -c x {x,y)\ < e(A) 
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and 

(4.9) 8(X)\c x (x,y) - &\(x, y)\ < e(A). 

We first show the bound (4.8). By Lemma 4.1 it is enough to show for 
all \x-y\ <5/\ l / d that <5(A)|E[£ A (z;:P A/t U y)£\(y,V\ K U x) - £\(x;V X k ( x ) U 
Uz)]| -*0 as A —> oo. 

We proceed as in the proof of (i), but now consider the event Cl xy that 
the radii of stabilization R(X^^x) and R(X 1 ^ d y) of £ for X l t d x and A 1 / d y, 
respectively, are both less than 5, that V\ K = Vx K (x) on the ball B^^ x y/d(x), 
and that V' Xk = V' Xk ^ on the ball B^j x y/d{y). Since y is within <5/A 1//rf of 
x, the probability that V Xk ^V\ K ( X ) on the ball B^ 5 ^ x y/d(y) is less than 
^d,8t K ((5/X) 1 ^) + ai5 ~ a2 . Therefore, under polynomial stabilization 
P[n c x J < 2u d St K ((8/X) 1 / d ) + 2oi<5-“ 2 . 

The triangle inequality, the moment condition (2.2) with p = 4 and Holder’s 
inequality give 

5|E[£ a (x; V Xk U y)£ A (y; V Xn U x) - £ A (x; V Xk ^ x) U y)£ A (y; V x ^ x) U x)} \ 

< a 3 <5(Ee 4 ) 1/4 P[^ J 3/4 < a 4 5(5t K «8/\) 1 / d ) + 6~ a *f\ 

We may similarly show 

5\c x (x,y) - c x (x,y) \ < a 3 8(E^) 1/4 P[n c x J 3/4 < a 4 5(St K ((S/X) 1/d ) + 8~ a *f 4 , 

which tends to zero as A —> oo since > 2. Thus (4.8) is satisfied. 

We now show the bound (4.9). Notice that the SV(^) assumption on £ 
implies automatically a somewhat stronger statement that the convergence 
is uniform not only on each fixed compact K , but also on the balls of radius 
S s ( A) —> oo as A —> oo. Even more strongly, we have convergence to zero with 
rate o((5 s (A))^ 1 ) uniformly on balls of radius <5 s (A)/A 1,/d , that is, 

sup (A) E[\^(y-V XK{x) ) - -0 

y: \x-y\<& s (\)/\ 1 / d 

as A —> oo. Thus by Holder’s inequality (p = 4, q = 4/3) we have 
(5 s (A)|E[£ A (x;P AK ( 3 .))£ A (y;P AK ( a .)) - €x(x;fi XK ( x ))€ x (y;P XK ( x ' ) )]\ 

< (A)(E[£ 4 (x-P X k{x) )} ) 1/4 (E|£ A (y; V Xk(x) ) - ££(y; V Xk{x) ) | 4/3 ) 3/4 - 0 
uniformly over balls around x of radius 8 S ( A)/A 1//rf . Similarly, 

5 a (A)|E[£ A (z; P Xk ( x ))&(v, V Xk {x)) ~ £aOg V X k(x))Z\(v, Pkw)]h0 
as A —> oo. Therefore the bound (4.9) holds on such balls. On the comple¬ 
ment of these balls, the bound (4.9) also holds by the polynomial decay of 
correlation functions given by Lemma 4.1. Therefore the bound (4.9) holds 
and the proof of Lemma 4.2 is complete. □ 
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4.3. Convergence of variance. We establish the convergence of A -1 Var [(/, fi XK )] 
for all / € C(A). Convexity and compactness of A implies the smoothness 
condition limn-^ n~ 1 d r (n 1 ^ d A) = 0 for all r > 0, where d r (n 1 / d A ) denotes 
the volume of the r-neighborhood of the boundary of n l ^ d A. Recalling (4.4), 
it suffices to show for all x € A distant at least 2(5/A) 1 ^ from dA , that for 
large A, 

(4.10) f(x)q x (x) + \f f(y)c x (x,y)n(y)dy 

Ja 

is close to 

(4.11) f(x)q x (x) + x[ f(x)c x x (x,y)n(x)dy. 

J 

Without loss of generality, assume Supp / is the set A. 

Lemma 4.2 (i) implies that for all x the difference of the first terms in 
(4.10) and (4.11) goes to zero as A —> oo. The difference of the integrals in 
(4.10) and (4.11) equals 

(4.12) A [ [c x (x,y)f(y)K(y)-c x x (x,y)f(x)K(x)]dy. 

J]SL d 

Let 6 := 5(A) be as in Lemma 4.2 and let B^ S / x y/d(x) be the ball of radius 

(< 5/A ) 1//d around x. To evaluate the integral (4.12), we integrate separately 
over B, s / x \i/d (x) and W l \ B (s / x y/d (x). The integral over B^ s / x y/d (x) involves 
the difference 

A / (c x (x,y)f(y)K(y)-c x x (x,y)f(x)K(x)]dy, 

which we split as 

A / (c x (x, y) - c x x (x, y))f(y)n(y) dy 

(4.13) + A / c x x (x,y)(f(y) -f(x))n(y)dy 

+ A [ c x x (x, y)f(x)(n(y) - k(x)) dy. 

•%/a) 1 /^ 

The hrst integral is bounded by the product of A, the volume of B^ s ^ x ^i/d(x) 
and the maximum of the integrand ( c x (x,y) — c x (x, y))f(y)n(y). However, 
since y is distant at least (5/A) 1 ^ from dA , the product goes to zero by 
Lemma 4.2(h). The second and third integrals also tend to zero as A —> oo 
by the bound (4.5) and Lemma 4.2 (iii). 

Since / and k are bounded, the integral in (4.12) over R. d \ B^i x y/d (x) is 
bounded by 


(4.14) 


cf [c x (x,y) + c x (x,y)\d(A 1/d y) 
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which by Lemma 4.1 is bounded by 

C [ (r(\z — x\/2)) 1 ^ 2 + (r(\z - x\/2)) 1 ^ 2 dz. 

jRd\B sl/d (x) 


The above integral is bounded by 2CuJd J^/d / 2 { rK {t)) 1 ^ 2 t d 1 dt which tends 
to zero as A —> oo by assumption. We conclude that (4.14) converges to zero 
uniformly for all x € A distant at least 2(5/X) 1 ^ from dA. Hence, 


. -l 


VarK/,/4)]- / n(x)f(x) 

J A 


f(x)q\(x) + A / f(x)c x x (x,x+y)n(x)dy 


dx 


converges to zero as A —> oo. The equivalences q\(x) = q(x) and c%(x, x + y) = 
c x (x, x + \ l t d y) yield (4.4) as desired: 


(4.15) 


A Var[(/,/4)J 
-> [ fix) 


f(x)q(x)+ / / (x)c x (x, y)n(x) dy 


k(x) dx. 


5. Proof of Theorem 2.4. We will only prove Theorems 2.4 and 2.5, 
since they are clearly a generalization of Theorems 2.1 and 2.2. We will first 
prove Theorem 2.4. We have already established Theorem 2.4(i) under the 
hypothesis that £ satisfies (2.2) for p = 4, and now to prove Theorem 2.4(h) 
we assume that £ satisfies (2.2) for all p > 0. 


5.1. Cumulant measures. Recall that C(A) denotes the continuous func¬ 
tions /: A —> R. To prove convergence of the finite-dimensional distributions 
of ,A > 1, in Theorem 2.4, it suffices to show for all test functions 

/ G C(A) that the Laplace transform of the random variable A” 1 / 2 (/, p,^ K ) 
converges as A —> oo to the Laplace transform of a normal random variable 
with mean zero and variance \ f A f 2 (x)V^(x,n(x))n(x)dx. In other words, 
it suffices to show for all / E C(A) that 

(5.1) lim Eexp(A -1//2 (—/,)) = exp \ [ f 2 (x)V^(x,n(x))n(x)dx . 

A—>oo I J A 

We will use the method of cumulants to show (5.1). We first recall the 
formal definition of cumulants. Let W := and formally expand (5.1) in a 
power series in / as follows: 

°° \~ 

(5.2) Eexp(X~ 1/2 (-f,fll K )) = l + J2 — 

k= 1 

where f k : R dk —► R ,k = 1,2,..., is given by f k (v i,. ..,v k ) = f(v i). f(v k ), 

and Vi € W, 1 < i < k. M k is a measure on M. dk , the kth moment measure 
(page 130 of [10]). 


k / 2 ((-f) k ,M k )) 

k\ 
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We have 

k 

(5.3) dM% = mx(vi,...,v k )Y[K(vi)d(\vi), 

i= 1 

where the Radon-Nikodym derivative ... ,Vk) is given by 


(5.4) 


m\{vv k ) :=E 




2—1 


and where given v±,... ,Vk we abbreviate notation and write for all 1 < i < 
k, £\(vi]V xk) for £\(vi]Vx K ) -E[£ x (vi-,V\ K )} and £(vf,Vx K ) for €(vnVx K U 
{vj}j =1 ). For each fixed k, the mixed moment on the right-hand side of 
(5.4) is finite uniformly in A by the moment bounds (2.2). Likewise, the kth 
summand in (5.2) is finite. 

When the series (5.2) is convergent, the logarithm of the Laplace func¬ 
tional gives 


log 


^ A~*/ 2 ((-/)*,M A fc ) 

h k '- 


^ a-*/ 2 ((-/M) . 

7! 

7=1 


the signed measures c l x are cumulant measures (semi-invariants [20] or Ursell 
functions). Regardless of the validity of (5.2), all cumulants c l x ,l = 1,2,..., 
admit the representation 

A= E (-ir'fp-iK'-E. 

Ti,...,T p 


where M x l denotes a copy of the moment measure on the product 

space W Ti and where Ti,... ,T p ranges over all unordered partitions of the 
set 1(see page 30 of [20]). More generally, c\ := c\(T) is the cumulant 
measure on W T with the representation 

cl= E (-1) p Ep-1)!^---mE 

Ti,...,T p 

where Ti ,..., T p ranges over all unordered partitions of the set T. The first 
cumulant measure coincides with the expectation measure and the second 
cumulant measure coincides with the covariance measure. The cumululants 
c l x ,l = 1,2,..., all exist under the moment condition (2.2). In what follows 
we make critical use of the standard fact that if the cumulants c l x of a random 
variable Z vanish for l > 3, then Z has a normal distribution. 

We will sometimes shorten notation and write M fc ,m and c l instead of 
M x ,m\ and c l x . 








GAUSSIAN LIMITS FOR RANDOM MEASURES 


31 


5.2. Cluster measures. Since c\ coincides with the expectation measure, 
we have ( f,c\) = 0 for all / G C(A). We already know from Section 4 
that A _1 (/ 2 ,c|) = A -1 Var[(/,/4j] -> f A f 2 (x)V^(x,K(x))K(x)dx. Thus, to 
prove (5.1), it will be enough to show for all k > 3 and all / G C{A) that 
as A —* oo (see, e.g., Lemma 3 of [28]). 

A cluster measure U?’ T on W s x W T for nonempty S,T C {1,2,...} is 
defined by 

U^' t (A x B) = Mf uT (A x B) - Mf (A)AfJ (B) 

for all Borel A and B in W s and W T , respectively. 

Let Si and S 2 be a partition of S and let T\ and T 2 be a partition of T. A 
product of a cluster measure U^ 1 ,Tl on W Sl x W Tl with products of moment 
measures on W S2 x W Ti will be called a ( S , T) semi-cluster measure. 

For each nontrivial partition (S, T) of {1,..., k} we next provide a repre¬ 
sentation of the /cth cumulant c k as 

c k = J2 a((S 1 ,T 1 ),(S 2 ,T 2 ))U s ^M^M^, 

(5i,Ti),(5 2 ,T 2 ) 

where the sum ranges over partitions of (1 consisting of pairings 

(S'ljTi), (. S 2 ,T 2 ), where S\ , S 2 C S and Ti,T 2 C T, and where a((S'i,Ti), (5 2 ,T 2 )) 
are integer-valued prefactors. 

In other words, for any nontrivial partition ( S,T) of {1,... ,k}, we show 
that c k is a linear combination of ( S,T ) semi-cluster measures. We were 
unable to find a proof of this in the literature and provide it. 

Lemma 5.1. For each nontrivial partition ( S,T ) of{l,...,k} we have 
ck = E a((S 1 ,T 1 ),(S 2 ,T 2 ))U Sl ’ Tl M^M^. 

(Si,Ti),(S 2 ,T 2 ) 

Proof. The proof involves some notation and definitions. The moment 
measures AP are expressed in terms of the cumulants via 

M j = E c(Ti)---c(T p ), 

Ti,...,T p 

where the sum is over all partitions of {1 ,,j}, that is, unordered collec¬ 
tions Ti, T 2 ,... ,T p of mutually disjoint subsets of {1,..., j} whose union is 
{1,..., j} (page 27 of [20], or [18]). Similarly, for any sets S and T, 

(5.5) M SuT = J2 c(S u Ti)---c(S p ,T p ), 

(SuBhiSp^p) 

where the sum is over all partitions of S' U T, where Si C S, Ti C T. A typical 
element (S*,Tj), 1 < i <p, of a partition thus involves a pair of sets, one a 
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subset of S and the other a subset of T. Some partitions of S U T are such 
that the empty set appears in each pair 1 < i <p. Call these the 

degenerate partitions. 

We now prove Lemma 5.1. Split (5.5) as 

(5.6) M Sut = J2c{- ■•)■■■£(■■■)+ 51 c(- ■■)■■■<■■■), 

{...} {■■■}* 

where {•••}* denotes degenerate partitions. The first sum contains the curnu- 
lant c(SUT) as well as products of lower-order cumulants, that is, cumulants 
of the form c(Si U Tj), where £* U Tj is a proper subset of {1,..., k}. Since 
each c(Si U Tj) is a product of moment measures, it follows that the first 
sum contains c(S U T) as well as linear combinations of ( S , T) semi-cluster 
measures. The second sum is just the product of M s and M T . Thus the 
cumulant measure c(S UT) is 

c(S U T) = M SuT - M s M t + l.c., 

where l.c. denotes a linear combination of (S, T ) semi-cluster measures. Since 
M SJT — M s M t is a (S, T) cluster measure, it follows that c(S U T) is a 
linear combination of semi-clusters. In particular, if ( S,T ) is a partition of 
{l,...,fc}, then c k is a linear combination of (S,T) semi-cluster measures. 
□ 


The following bound is critical for showing that A c k ) —> 0 for k > 3 

as A —> oo. 

Lemma 5.2. If £ is exponentially stabilizing, then the functions m\ clus¬ 
ter exponentially; that is, for positive constants Ajj and Cfj one has uni¬ 
formly 

|m A (xi, . . m x (x i, . . .,Xj)mx(yi, ■■■,yi)\ 

< A jt iex.p(-Cjj5\ 1/d ), 

where 6 := mini<j< y .i< p </1 Xi — y p \ is the separation between the sets (xi) 3 i=l 
and (y p ) l p=1 . 

Proof. The proof is similar to that of Lemma 4.1. With probability at 
least 1 — exp(— dX 1 ^/C), the radius of stabilization for each A l ^ d Xi, 1 <i <j, 
and each \ l / d y p , 1 <p <1, is less than A 1//d <5. Let Ejj := Ejj(d) denote the 
event for which all such radii are less than A 1 ^^. On E,jj the stabilization 
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balls do not intersect and therefore 


E 


j 

[ 

.i =1 


W_^\{xi\V XK ) n €\(y P ','P\K)iE j ' l 

p =1 


= E 


3 

[ 

i =1 


~l^x(Xi\Vx 


E 


II €\(y P \'P\K)i-E j , l 

p =1 


Holder’s inequality and the moment conditions imply there is a constant 
Ajj such that 

\m x {xi, ...,x j ,y 1 ,...,yi)~ m x (x i,... ,Xj)m x (yi, . ..,yi) \< A jti {P[FF u ]) 1/2 . 
Since decays exponentially in <5, Lemma 5.2 follows. □ 


The next lemma specifies decay rates for the cumulant measures. Such 
decay rates are useful in establishing moderate deviation principles and laws 
of the iterated logarithm for the measures Ji\ h . [4]. Here we simply use the 
decay rates to conclude the proof of Theorem 2.4. 


Lemma 5.3. For all f £ C(A) and for all k = 2,3, ..., we have A k ^ 2 {f k ,c k ) 

0(ll/llLA (2 - fc)/2 ). 


Proof. We need to estimate 

x -k /2 I f( Vl )... f( Vk )dc x (v 1 ,...,v k ). 

J A k 

Let A*, denote the diagonal in W k , that is, v\ = V 2 = • ■ • = v k . For all v := 
(v\ ,... ,v k ) € A k , let D k (v) denote the distance to the diagonal. 

Let n(fc) be all partitions of {1,2,..., k} into two subsets S and T. For all 
such partitions consider the subset u(S, T ) of A s x A T having the property 
that v € a(S,T) implies d(x(v),y(v)) > D k (v)/k, where x(v) = v D A s and 
y(y) = »(1 A T . Since for every v := (vi,... ,v k ) £ A k , there is a splitting 
x := x(v) and y := y(v) of v such that d(x,y) > D k (v)/k, it follows that A k 
is the union of the sets cr(S, T ), (S', T) £ n(fc). The key to the proof of Lemma 
5.3 is to evaluate the cumulant c\ over each a(S,T). We then use Lemma 
5.1 and adjust our choice of semi-clusters there to the particular choice of 
(S,T). 

By Lemma 5.1, the cumulant measure dc k (v i,..., v k ) on cr(S, T ) is a linear 
combination of (S, T) semi-cluster measures of the form 

£ a((Si,ri),(S 2 ,r 2 ))t/‘ Sl ’ Tl M |52| M |T21 , 

(5i,T 1 ),(5 2 ,T 2 ) 
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where the sum ranges over all partitions of {1,..., k} consisting of pairings 
(Si,Ti), (S 2 ,T 2 ), where S\,S 2 C S and Ti,T 2 C T, and where a((5i,Ti), (S 2 ,T 2 )) 
are integer-valued prefactors. 

Let x and y denote elements of A 5 and A T , respectively. Let x and y 
denote elements of Tl 51 and A Tl , respectively, and let x c denote the comple¬ 
ment of x with respect to x and likewise with y c . The integral of / against 
an (S, T) semi-cluster measure has the form 

A “ fc/2 [ f(v 1 )---f(v k )d(M[ S 2 \x c )Ui +j (x,y)Mx 2 l (y C ))- 

Ja{S,T) 

Letting u\{x,y) := m x (x,y) — m x (x)m x (y), and recalling (5.3), the above is 
bounded by 

k 

(5.7) X~ k/2 f f(v k )m x {x c )u\(x,y)m\{y c ) k(u;) d(Au;). 

J*{S,T) f=i 

Decompose the product measure Yli=i K ( v i) d(Xvi) into two measures, one 
supported by the diagonal A k and the other not. Off the diagonal, the 
integral (5.7) is bounded by 

D\\f\\ k oo^~ k/2 / exp (-Ct)P[D k >t]dt = 0( A" fc / 2 A), 

Jo 

since u\ decays exponentially with the distance to the diagonal (Lemma 
5.2), the mixed moments m\ are uniformly bounded, and since the differ¬ 
ential of a volume element of points at a distance greater than t from the 
diagonal is bounded by the Lebesgue measure of the diagonal. Integrating 
over the diagonal measure Xk(v\) dv\, and using the boundedness of /, we 
thus bound (5.7) by D||/||^ 0 A _i; / 2 A for some constant D. Since this esti¬ 
mate holds for all a(S, T), (S, T) e Ll(fc), where A k is the finite union of sets 
a(S,T), Lemma 5.3 holds. □ 

6. De-Poissonization: proof of Theorem 2.5. De-Poissonization involves 
a significant modification of de-Poissonization arguments for CLTs for translation- 
invariant functionals (Section 4 of [24]) defined over homogeneous point sets 
and thus we provide the details. 

Before de-Poissonizing, we need the following definition. For any X C A 
and / €E C(A), let 

(6.1) i#(*):=£/(z)M®;*)- 

xex 

Letting X m be a point process consisting of m i.i.d. random variables with 
density k on A, set R m}n := Hf(X m+1 )-Hf(X m ). 

The following coupling lemma is inspired by and follows closely Lemma 
4.2 in [24], 
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Lemma 6.1. Suppose f is exponentially stabilizing for k, and suppose H 
is strongly stabilizing. Let e > 0. Then there exists d > 0 and no > 1 such that 
for all n>no and all m,m' £ [(1 — 5)n, (1 + <5)n] with m < m', there exist 
random variables X,X' with density k and a coupled family of variables 
D:=D(X), D' :=D(X'),R:=R(X,X'),R' :=R'(X,X') with the following 
properties: 

(i) D and D' each have the same distribution as f(X) A^(X, k(X)); 

(ii) D and D' are independent; 

(iii) (R,R r ) have the same joint distribution as (R m ,n, Rm 1 ,n)! 

(iv) P[{\D - R\ > e} U {|Df - R !| > e}} < e. 

Proof. We will modify the proof of Lemma 4.2 of [24], Suppose we 
are given n. Let X^X',Y\,Y 2 ,... be i.i.d. random variables with density k 
on A. On the probability space (f l,P,P), let V := V nK and V' := V' nK be 
independent Poisson processes on A with intensity measure nn{x) dx. 

Let V" be the point process consisting of those points of V which lie closer 
to X than to X' (in the Euclidean norm), together with those points of V 
which lie closer to X' than to X. Clearly V" is a Poisson process also having 
intensity measure nn(x) dx on A and, moreover, it is independent of X and 
ofX’. 

Let N denote the number of points of V” (a Poisson variable with mean 
n-volA). Choose an ordering on the points of V", uniformly at random from 
all N\ possible such orderings. Use this ordering to list the points of V" as 
Wi, W 2 ,..., W N . Also, set Wn+i = Y \, Wv +2 = Y 2 , W N+3 = Y 3 and so on. 
Let 

R := R(X,X 0 := H^Wt ,..., W m ,X}) - H^Wu ..., W m }) 

and 

R' ■= R'(X,X') X, X'}) - Hf({W u ..., W m ,_ u X}). 

X, X', Wi, W 2 , W 3 ,... are i.i.d. variables on A with density k, and there¬ 
fore the pairs (R,R') and (Rm,n, Rm 1 ,n) have the same joint distribution as 
claimed. 

For all x £ and r > 0, let B(x,t) denote a stabilization ball for £ at x 
with respect to V T . Recalling A X (A?) := H(X Lli) — H(X) we put 

D(x) := f(x) A x ('P k(x) n B(x, k(x))) 

and 

D '(x') := f(x') A x fV' K{xl) n B{x\ n(x'))). 

Let D := D(X), D’ := D'(X'). Then D and D’ are independent, and by 
strong stabilization of H , given (X,X') = (x,x'), have the same distribution 
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as f(x)A^(x, k(x)) and f(x') A*(x', k(x')), respectively. It remains to show 
that P[{\D - R\ > e} U {\D' - R'\ > e}] < e. 

Without loss of generality, we may couple V K t x \ and n x l d V nK such that 
for all Borel sets Bed 

P[Pn K (B) ^ V nK ( x )(B)\ <n f \n(y) - n(x)\dy. 

J B 

Given (X,X') = (x,x r ), for every K > 0 and n= 1,2,..., let Fx. n '■= 
F K , n (x) be the event that B{n}/ d x, n(x)) C Bx{n L ^ d x') and that V K ^ and 
n)/ d VnK coincide on the ball Bxip^^x). As in the proof of the limit (2.3), 
we have by the uniform continuity of k and the Lebesgue point property of 
x that for K and n large enough, P[F^ ] < e/9 uniformly in x. Similarly, 

let Fx, n ■= P'Kn( x ') be the event that B(n l / d x', k(x')) C Bx(n ] / d x') and 
that V K ( x i) and r^^ d V' nK coincide on Bxi^^x'). For K and n large enough 
P[(F' Kn ) c } < e/9 uniformly in x'. 

Thus given (X, X') = (x,x'), on sets Fx, n and F k n of probability at least 
1 — e/9, we have 

D(x) = f{x)A nl/dx {n 1,d V nK n B(n 1/d x , k(x))) 

and 

D'(x') = f(x')A nl/dx ,(n 1/d V' nK n B{n l/d x', k(x'))). 

Thus given (X,X') = (x,x'), we need only show that 

P[\f(x)A n i /dx (n 1/(l V nK nB(n 1/d x,n(x))) - R(x,x')\ > e] < e 

and 

P[\f(x')A nl / dx ,(n 1/d V nK nB(n 1/d x',K(x'))) - R 1 {x,x')\ > e] < e. 

We will show the first bound only; the proof of the second bound is identical. 

We now follow [24], page 1018. Choose K large enough such that P[S > 
K] < e/9. For all w G and r > 0, let Q r (w) = [—r,r] d + w be the cube 
centered at w. Let tf denote the modulus of continuity of / € C{A) and 
find b(n) such that 6(n)t/(iv/n 1//d ) —>• 0. Given e and K as above, let n> 
n(e, K) be so large that 

(6.2) supb(n) ■t f (K/n 1/d )E\£(n 1/d W 1 ;n 1 / d P nK Un 1/d x)\ < e 2 /72 

X 

and 

(6.3) b(n) ■ tf(K/n 1 / d )E\l;(n 1 / d Wi]n 1 / d V nK )\ < e 2 /72. 

Take n so large that except on an event (denoted Eq) of probability less 
than e/9, the positions of ii}/ d x and ii}/ d x' are sufficiently far from the 
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boundary of n l ^ d A and from each other, that the cubes Qxi^^x) and 
Q K (n 1/d x') are contained entirely within v}^ d A (possible by the regularity 
of dA), and also are such that every point of Qk^^x) lies closer to v}^ d x 
than to v}/ d x' and every point of Qk(j i^^x') lies closer to v}^ d x' than to 
•n}/ d x. 

Set 5 := e(2/i) _d /18. We assume \m — n\< 5n and | nn! — n\ < 6n. For n 
large enough, except on an event (denoted E\) of probability at most e/9, 
we have \N — m\ < 25n = e(2it) _d n/9, and likewise \N — m'\ < e(2it) _d n/9. 

Let E be the event that the set of points of n 1//rf {kFi,..., W m } lying in 
Qk^^x) is not the same as the set of points of V lying in Quin^^x). 
This will happen either if one or more of the (N — m) + “discarded” points 
of n}-l d V" or one or more of the (m — N) + “added” points of n 1 ^ d {Yi,Y 2 ,...} 
lies in Qxin^^x). For each added or discarded point, the probability of lying 
in (5if(n 1//rf x) is at most (2 K) d /n, and so the probability of E, given that 
E\ does not occur, is less than e/9. 

We now compute 

P[\D(x) — R(x,x')\ > e] 

<P[\D(x) — R(x,x') |lF Xjn nugnEfnE c n{S<A} > e / 2 ] 

+ P[\D(x) - -R(®,T / )|l(F Jf , n nEgnEfnu c n{S<A}) c > e/ 2 ] 

<P[\D(x) — R(x,x')\l FKnnE c nE c nE e V { S<K y >e/2] 

+ P[Fk,n] + P[E o] + P[Ei] + P[E \ El] + P[S > K]. 

The last hve terms are bounded by e/9 for large n. Now consider the 
first probability. On the set Fk,u 0 Eq 0 Ef 0 E c n {S < K} the difference 
| D(x) — R(x,x')\ equals 

\f{x)A n i/ dx (n 1/d 'P nK 0 B(n 1/d x, k(x))) - R(x, x')\, 
which by strong stabilization of H is bounded by (since S < K ) 

< E l/(Wi)-/(x)||^(n 1/d Wi;n 1 / d P rM( Un 1 / <i s)| 

n 1 / d W i eQ K (n 1 / d x) 

+ E I HWi) - fixmin^W^n^VnJl 

n 1 / d Wi£Q k ( nl I d x ) 

By definition of tf, the above is bounded by 

<t f (K/nV d ) E |e(n 1/d VF i ;n 1 / d iP nK Un 1 / ,i x)| 

n}/ d Wi £Qk 

+ t f {K/n 1 / d ) E le(n 1/d W i ;n 1 / d iP nK )|. 

n 1 / d WiGQfc(n 1 / <i ai) 
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Let Nk := card(n 1/,rf {VLj} n Qk^^x)). Then the first term in the above is 
bounded by 

t f (K/n 1/d ) J2 \^(n 1/d Wi;n 1/d V nK Un 1/d x)\l NK < b{n) 

71 1 / d Wi€Q Kip* 1 / d x) 

+ t f (K/n 1/d ) J2 \^(ri 1/d Wi-,n 1/d V nK Un 1/d x)\l NK>b{n) , 

n 1 / d Wi£Qfc(n 1 / d x) 


with a similar bound for the second term. Therefore, combining all of the 
above bounds 


P[\D(x) — R(x,x')\ > e] 

< 5e/9 + P[\D(x) - R(x,x’)\l FKnnE c nE c nE c n ^ s<K > > e/2] 

b{n) 


(6.4) < 5e/9 + P 


tf(K/n l l d ) ^2 \€(i r i' 1 / d Wi;n 1 / d V nK U n 1 ^ d x)\ > e/8 


+ P 


b(n) 


i —1 


tfiK/n 1 ^) J2 \C(n 1,d Wi-n l,d V nK )\ > e/8 


i= 1 


+ 2 P[N>b(n)]. 


Using Chebyshev and the bounds (6.2) and (6.3), the second and third terms 
in (6.4) are bounded by e/9 for n large enough. For n large, the last term 
in (6.4) is bounded by e/9, since N is a.s. finite. 

Now integrate over all pairs (x,x f ) to obtain the desired result. □ 


The next lemma extends Lemma 4.3 of [24], 

Lemma 6.2. Suppose £ is exponentially stabilizing and satisfies (2.2) 
for all p > 0. Suppose H is strongly stabilizing and satisfies the bounded 
moments condition for n. Let (h(n)) n > i be a sequence with h(n)/n—* 0 as 
n —> oo. Then 

(6.5) lim sup \ER m , n -¥.[f(X)At(X,K{X))]\ = 0. 

n ~^°° n—h{n)<m<n+h{n) 

Also 

(6.6) lim sup \ER mtn R m ,, n -E[f(X)AZ(X,K(X))] 2 \ = 0 

n— ^°° n—h(n)<rn<m'<n-\-h(n) 

and 

(6.7) lim sup < oo. 

n ~n—h(n)<m<n+h(n) 

Proof. We will follow the proof of Lemma 4.3 of [24], Let m be an 
arbitrary integer satisfying n — h{n ) <m<n + h(n). Let e > 0. Provided 
n is large enough, by Lemma 6.1 we can find coupled variables D and R, 






GAUSSIAN LIMITS FOR RANDOM MEASURES 


39 


with D having the same distribution as f(X)A^(X, k(X)), with R having 
the same distribution as R m ,n, and with P[\D — R\ > e] < e. Then 

E-R m , n = E. R = E[L>] + E [(R — D)\\ R _ D \ >e ] + E \(R — D)l\ R _ D \< e ]. 

By Cauchy-Schwarz, the moments condition (2.2), and the fact that the 
bounded moment condition implies E[D 2 ] < oo (Lemma 4.1 of [24]), we 
have 

n{R-D)l\ D _ R>£ \<Ce 1 ' 2 . 

Since e is arbitrarily small, (6.5) follows. The proof of (6.7) is similar and is 
omitted. 

Next we consider m, m! with n — h(n ) <m <m! < n + h(n). By Lemma 
6.1, there are coupled variables D,D',R,R' such that D and D' are inde¬ 
pendent and each has the same distribution as f(X)A^(X, k(X)), (R,R r ) 
have the same joint distribution as (Rm,n, Rm',n ), and 

P[{\D - R\ > e} U {| D' - R'\ > e}} < e. 

AowElRR'j-ElDD 1 ] =E[R(R! - D')]+E[D’(R-D)]. By Cauchy-Schwarz, 
we again obtain the bounds E[i? | R' — D' |] < Ce 1 / 2 and K[D'\R — D\] <Ce 1//2 . 
It follows that the difference E[i2i? ; ] — E [DD'\ can be made arbitrarily small 
and (6.6) follows. □ 


Proof of Theorem 2.5. We first prove the limit (2.18). Given / £ 
C(A), X C. A and recalling (6.1), let H } { := H^(X n ) and H'{ := H[(T nK ). 
Assume that V nK is coupled to X n by setting V nK := {X\,X 2 ,... ,X]\r n }, with 
N n an independent Poisson variable with mean n. 

To prove (2.18), it is enough to show for all / £ C(A) that 


( 6 . 8 ) 


H f - E H f 

_n 

n 1 / 2 


N(0,rj), 


where N(0,r 2 ) denotes a mean zero normal random variable with variance 


(6.9) 

Letting a := 

( 6 . 10 ) 


t 2 := I f 2 (x)V^(x,K(x))K(x)dx 


IA 


f(x) E[A^(x, k(x))]k(x) dx 


E[/(A)A^(A, k(A))], the first step is to prove that as n —> oo, 
E -H^ {.N n - n)a) 2 } 0. 


To do this, we employ the coupling Lemma 6.1 and follow pages 1019 and 
1020 of [24] verbatim. The second step is to prove that 


( 6 . 11 ) 


, Var Hf 
lim - - 


n 
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However, this follows from the identity 

nT x l 2 H'£ = + n _1 / 2 (jv n — ri)a + n~ 1,/2 (f/X - Hl — ( N n — n)a). 


The third term in the above has variance tending to zero by (6.10); the 
second term has variance a 2 and is independent of the first term. Letting 
aj := J A f 2 (x)V^(x,K(x))n(x)dx, it follows that 


o 


= lim 


Var [H\ 


'J) 


n 


= lim 

71—XX) 


Var (Hf) 


a 


n 


that is, (6.11) holds. The limit (6.8) follows as on page 1020 of [24], thus 
establishing (2.18). 

Let us now show (2.19). The above shows that the sequence of distribu¬ 
tions rt _1 / 2 (p^, /) tends to a limiting normal random variable N(0,rj) for 
every / £ C(A). Taking f = fi + f 2 and using simple algebra shows that 
the limiting Gaussian field has the desired covariance matrix (2.19). This 
completes parts (i) and (ii) of Theorem 2.5. 

To prove Theorem 2.5(iii), it suffices to show that 


( 6 . 12 ) 


lim Var [H n {X n )\ 

n —xx) 77, 


= [ V^x,K(x))K(x)dx 
Ja 


j E[A^(x, k(x))]k(x) 



> 0 . 


We accomplish this by modifying the approach in Section 5 of [24]. 

We write H n (X n ) — E H n (X n ) as a sum of martingale differences as fol¬ 
lows. Let = a(Xi :n ,... and write Ej for conditional expectation 

given Xi. Define martingale differences D{ := E iH n (X n ) — Ej_iiL„(fL n ). Then 
H n (X n ) - E H n (X n ) = E”=, Di and 


V a r[H n (X n )] = J2nDf}. 

i= 1 

It suffices to show that there exists a constant G > 0 such that for all 1 < 
i < n, E [Df] > C. 

Given i A nr, let — H n i^X rn ) H n (^X m \ { X$y ). Let G^ m — H n (^X m ^.\ \ 

{Xi}) — H n (X m \{Xi}). Then D { = E^G^-G^]. We set a : = E[A*(X, «(*))] 
and using nondegeneracy, take <5 > 0 such that P[A^(X, k(X)) > a + 4(5] > 
45. 

Define /: R —> R by f(x ) = 0 for x < a + 5 and f(x) = 1 for x > a + 25, 
interpolating linearly between a + 5 and a + 25. Let Y) := /(Ej[Gj jn ]). The 
remainder of the proof consists in showing that for n large and for i > 

(1 — £ 3 ) 71 , we have 

E[(Gj >n — a)Yi] > 45 2 and E[(G iiTl - a)V] < 25 2 . 
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These inequalities follow from Lemma 6.1 and pages 1021 and 1022 of [24]. 
Taken together, this implies for large n and i > (1 — £3 )n, that E[(Gj in — 
Gi, n )Yi] > 2 5 2 . Using the fact that 1/ is .^-measurable and lies in the range 
[ 0 , 1 ], we obtain 

26 2 <E[Y i E i (G ijn -G i , n )\<E[\E i {G i , n - G i>n )\] =E[\D i \], 

and hence, E [Df] > [E|Dj |] 2 > 4 5 4 > 0. Thus (6.12) holds, completing the 
proof of Theorem 2.5. □ 
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